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ABSTRACT 



The latest observations of line and continuum spectra emitted from the extended narrow line region (ENLR) of the Seyfert 2 galaxy 
NGC 7212 are analysed using models accounting for photoionization from the active nucleus and shocks. The results show that rela- 
tively high (500-800 km s~^) shock velocities appear on the edge of the cone and outside of it. The model-inferred AGN flux, which is 
lower than 10~^^ photons cm~^ s~^ eV~^ at the Lyman limit, is more typical of low-luminosity AGN, and less so for Seyfert 2 galaxies. 
The preshock densities are characteristic of the ENLR and range between 80-150 cm~^. Nitrogen and sulphur are found depleted 
by a factor lower than 2, particularly at the eastern edge. Oxygen is depleted at several locations. The Fe/H ratio is approximately 
solar, whereas the Ne/H relative abundance is unusually high, 1.5-2 times the solar value. Modelling the continuum spectral energy 
distribution (SED), we have found radio synchrotron radiation generated by the Fermi mechanism at the shock front, whereas the 
X-rays are produced by the bremsstrahlung from a relatively high temperature plasma. 
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1. Introduction 



Morphological and spectroscopic studies of Seyfert galaxies 
have shed light on the nature of active galactic nuclei (AGN). 
Specifcally, the spectra of type 2 Seyferts are generally domi- 
nated by strong forbidden and permitted lines revealing the phys- 
ical properties of the extended narrow line region (ENLR). The 
spectra indicate that radiation from the AGN dominates the radi- 
ation field; however, line ratios corresponding to low ionization 
levels and neutral lines could be explained by underlying shock 
wave hydrodynamics. The full width at half maximum (FWHM) 
of the lines reveal velocities up to 1000 km s"^ . Modelling of line 
and continuum spectra observed in many single regions through- 
out the ENLR provided a new dimension to our knowledge. Eor 
instance, the complex nature of NGC 7130, in terms of the in- 
termingled activity of starbursts, shocks, and of an active nu- 
cleus and their mutual location, could be traced by Contini et 
al. (12002a"). In such objects, e.g. NGC 4388 (Ciroi et al. 2003) 
or Mrk 298 (Radovich et al. 2005), the relative importance of 
the radiation flux from the AGN, of shocks, and of radiation 
from starbursts could be determined in the difl'erent regions of 
the galaxy. This method, which is useful for isolated galaxies, 
becomes fundamental for the analysis of Seyfert nuclei in merg- 
ing systems originating from collisions, even in case of multiple 
(generally double) nuclei that are not yet clearly identified as 
such by the observations. Eor instance, the biconical structure of 
the high excitation region which emerges from the toroidal ob- 
scuring material, can be fragmented by collision, and patches of 
highly ionized material can appear at the edges of the ENLR, 
as e.g. in NGC 3393 (Cooke et al. 2000). The interaction of the 
ENLR galactic matter with fast shocks, star formation, dust grain 
destruction etc. are amply analysed and summarised by Ramos 



Almeida et al. (l2009l) for NGC 7212 and other Seyfert 2 galax- 
ies. 

NGC 7212 (z=0.0266) belongs to a compact group of in- 
teracting galaxies, with two galaxies in a clear ongoing merger 
(Wasilewski 1981). Spectropolarimetry studies show the pres- 
ence of a hidden broad line region (BLR) (Tran & Kay 1992), 
clumpy nuclear morphology, and irregular diff'use emission. 
Munoz Marin et al. (12007b also claim that these galaxies man- 
ifest some cases in which in addition to the extended emission, 
the UV light stems from the knots or clumps likely produced 
by star clusters. They identify the bright clumps to the south as 
star-forming regions with a certain contribution from the AGN. 
Exploring whether photoionization or shocks dominate in NGC 
7212, Bianchi et al. (120061) found a striking resemblance of the 
[O III] structure with soft X-ray emission, a clear indication that 
shocks are at work. 



In an accompanying paper, Cracco et al. (12011b presented 
new observations of NGC 7212 and suggested that its ENLR 
gas has an external origin, likely due to gravitational interaction 
eff'ects in act in this triple system. The optical spectra detailed 
by Cracco et al. are rich in number of lines: e.g. oxygen from 
three ionization levels ([OIII]5007, [OII]3727 and [OI]6300), 
He II 4686 and He 1 5876, [Ne III]3869, [N II]6548,6583, the two 
[SII] lines 6716 and 6731, lines from higher ionization levels, 
e.g. [ArIV]4713 and [EeVII]6087, in addition to 11/3 and Ua. 
The line ratios constrain the physical conditions of the emitting 
gas and the relative abundances. These authors pointed out for 
the first time the presence of an extended ionization cone in NGC 
7212, with high values of [O III]/Hyg, up to 3.6 kpc. In fact, the 
Veilleux & Osterbrock (11987b diagnostic diagrams indicate that 
the main source of ionization, at least in the entire field of view 
of the integral field spectra, is the AGN. However, high values of 



1 



Contini et al.: Depletion of heavy elements in NGC 7212 



[N II]/HQf and [S IIJ/Hof are observed towards the region of grav- 
itational interaction between the interacting galaxies, suggest- 
ing a possible combination of ionization by the active nucleus 
and shocks. On the other hand, no contribution from stellar ion- 
ization is required to explain the observed emission line ratios. 
An archival broad-band Hubble Space Telescope (HST) image 
of NGC 7212 obtained with the F606W filter revelas a struc- 
tur e com posed of clouds or filaments (see fig. 19 in Cracco et 
al. 1201 li) . The evidence of dust located in the ENLR also sup- 
ports the idea that the gas in the cone has a filamentary struc- 
ture. In addition, comparing the more significant line ratios (e.g. 
[OIII]/HyS, [OII]/HyS, [NII]/Ha, etc) with diagnostic diagrams, 
Cracco et al. found depletion of heavy elements such as N and 
O relative to solar. 

The emission line profiles of [OIII] inside the ionization 
cone display blue wings in the northern side of the cone, and 
red wings in the southern side, while in the nucleus the profiles 
are symmetric. This suggests the presence of gas in radial mo- 
tions, which is confirmed also by the analysis of high resolu- 
tion echelle spectra characterized by multiple kinematical com- 
ponents at diff'erent velocities. If radial motions are in act in the 
ENLR of NGC 7212, shock ionization is expected to play a non 
negligible role. 

Moreover, it is well known that the higher the flux from the 
photoionizing source, the stronger the [OIII] lines are, while 
emission lines such as [Oil], [Nil], etc., namely low ioniza- 
tion level lines, are significant when the gas is aff'ected by strong 
shock. Consequently, we decided to reproduce the observational 
data presented by Cracco et al. (2011) by means of models ac- 
counting for the coupled eff'ect of photoionization and shock, im- 
proving their interpretation, and adding information concerning 
the physical and chemical conditions throughout the galaxy. 

The structure of the paper is as follows: observations of the 
optical spectra are described in Sect.[2l In Sect.[3]the models and 
the modelling method are explained. Modelling results of the 
line and continuum spectra are presented in Sect.|3]and in Sect.[5l 
respectively. Discussion and concluding remarks follow in Sect. 
[6] and Sect.|7l 



2. Observations 

Cracco et al. (|2011|) present integral field spectra taken with the 
Multi Pupil Fiber Spectrograph (MPFS; Afanasiev et al. 2001 ) 
at the Special Astrophysical Observatory (SAO RAS). 256 si- 
multaneous spectra arranged in an array of 16 x 16 spatial ele- 
ments were obtained in the range 3800-7300 A, with a resolu- 
tion of 6 A. The field of view was 16x16 square arcsec, seeing 
was about 1.5''. At the distance of the galaxy the spatial scale is 
0.513 kpc arcsec"^ (Hq = 75 km s"^ Mpc"^), yielding field of 
view of 8.2x 8.2kpc^. 

The data were reduced with p3d (Sandin et al. 2010) and cor- 
rected for reddening, atmospheric absorption bands, and atmo- 
spheric refraction with iraiQ. The stellar component was sub- 
tracted with STARLIGHT (Cid Fcmandcs et al. 2005, 2007) and 
finally the spectra were measured with pan (Peak Analysis). 
Spectra are numbered from 1 to 256 and their positions in the 
field of view are reported in Fig.[TJ For further details about re- 
duction and calibration procedures, see Cracco et al. (1 2011b . 



^ iRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Fig. 1. 2D map of the stellar continuum at 5500 A. The numbers 
identify the spectra in Fig.fTJ The image size is 16 x 16 square 
arcsec. North is up and East is to the left. 



In Fig. [2] we show a few significant observed line ratios 
throughout the observed region in order to obtain a first hint 
about the ionization conditions of the emitting gas. The diagrams 
of the left column describe the spatial distribution of [OIII], 
[O II], and [O I] line ratios to HyS. It can be noticed that photoion- 
ization from the AGN dominates in the central regions where 
the [O III]/HyS and the [0 11]/UJ3 show smooth opposite trends. 
However, throughout all the ENLR, [OI]6300/H/5 fine ratios 
show moderate fluctuations (Fig. [2]) ranging between ~ 0.3-0.6. 
These line ratios are relatively low, implying that the spectra cor- 
respond to matter bound clumps. Such a narrow range is unusual 
in matter fragmented by the underlying turbulence which accom- 
panies relatively high velocity shocks. Eventually, [0 1]/UJ3 can 
be about 0, where the emitting cloud is so thin that there is not 
enough gas at temperatures < 10^ K. The [O I]/H/5 line ratios ob- 
served in NGC 7212 may suggest that the contribution from the 
galactic interstellar medium (ISM), that is ionized gas not per- 
turbed by high velocity shocks, emerges at some locations in the 
ENLR. 

In this work, we have tried to fit the [O I] lines consistently 
with [O III] and [O II] . Only in a few cases did the fit of the three 
oxygen ionization level lines to HjS require certain ISM contribu- 
tion, even when accounting for the typical observational errors 
of 0.12, 0.18 and 0.32, for [OII]/HyS, [OIII]/HyS and [OI]/HyS, 
respectively. 

3. The models 

Complex line profiles, as those observed in each position of 
NGC 7212, indicate a complex network of emitting clouds. The 
models help to disentangle the contributions of the diff'erent heat- 
ing and ionizing sources (e.g. photoionization from the AGN, 
diff'use secondary radiation, shock waves, etc.) to each spectrum 
by reproducing a large number of observed line intensity ratios. 
In this way we obtain the physical conditions prevailing in each 
position. 
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Fig. 2. Maps of the line ratios to HyS. From left to right: [O III] , [Fe VII] and [N II] (top) ; [O II] , [S II] 67 1 6 and He II (middle) ; [O I] , 
[S II]6731 and He I (bottom). Spatial scale is 1 arcsec px"^. North is up and East is to the left. 



3.1. Analysis of the spectra 

The clumpy and irregular morphology of NGC 7212 and the rel- 
ative large FWHM of the spectral lines indicate that in addition 
to the photoionization flux from the AGN, shocks are heating 
and ionizing the ENLR gas. We choose the shock velocity on 
the basis of both a close approximation to the observed FWHM 
of the line profiles and the best fit to the observed line intensi- 
ties. Actually, we consider that the FWHM represents roughly 
the shock velocity. So, when diff'erent profiles are intermingled 
in one spectrum and show a broad component compared with 
the bulk of the narrow profiles observed throughout the ENLR, 
we present both the spectra calculated by high and low velocity 
shocks. 

The physical parameters are combined throughout the cal- 
culation of forbidden and permitted lines (see Osterbrock 1989) 
emitted from a shocked nebula. 

In pure photoionization models, the density n is constant 
throughout the nebula, while in a shock wave regime, n is cal- 
culated downstream by the compression equation in each single 
slab. Compression depends on n, the magnetic field B and the 
shock velocity Vs . In models accounting for the shocks, both 
the electron temperature Tq and the electron density Hq are far 
from constant throughout each cloud, showing the characteristic 
profiles in the downstream region (see Sect. l4.2.5l) . So, even so- 



phisticated calculations reproduce approximately the highly in- 
homogeneous conditions of the gas, giving rise to some discrep- 
ancies between the calculated and observed line ratios. 

The ranges of the physical conditions in the nebula are in- 
ferred, as an initial guess, from the observed lines (e.g. the shock 
velocity from the FWHM) and from the characteristic line ratios 
(e.g. Hq and Tq). The physical parameters revealed by the obser- 
vations indicate that a steady state situation can be applied (Cox 
1972) in NGC 7212. In this case, the time t required for a par- 
cel of gas to cross the cooling region from the shock front to 
the recombination zone, for shock waves with = 100 kms~\ 
is found to be approximately 1000/^e yr (calculated by the re- 
combination coefl&cients), such that for an electron density Hq = 
100 cm"^, t = 10 yr. Shock velocities are not likely to change 
appreciably in such a short time, and therefore the steady state 
calculation may be regarded as adequate. 



3.2. The code: input parameters 

In this paper, the line and continuum spectra emitted by the 
gas downstream of the shock front are calculated by suma (see 
http://wise-obs.tau.ac.il/~marcel/suma/index.htm for a detailed 
description). The code simulates the physical conditions in an 
emitting gaseous cloud under the coupled eff'ect of photoion- 
ization from an external radiation source and shocks. The line 
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and continuum emission from the gas are calculated consistently 
with dust-reprocessed radiation in a plane-parallel geometry. 

The input parameters which refer to the shock are: the shock 
velocity Vs , the pre-shock density no, the pre-shock magnetic 
field ^0- The input parameters related to the radiation field are 
the power-law flux from the active nucleus F in number of pho- 
tons cm"^ s"^ eV"^ at the Lyman limit, if the photoionization 
source is an active nucleus, and the spectral indices ajjy = -1.5 
and ax = -0.7. F is combined with the ionization parameter 
U by U = (F/(n c (a - 1)) ((^h)"^^^ - (^c)"^^^) (Contini & 
Aldrovandi 1983), where Eu is H ionization potential and Eq is 
the high energy cutofl', n the density, a the spectral index, and c 
the speed of light. The secondary difl'use radiation emitted from 
the slabs of gas heated by the shocks is also calculated. The flux 
from the active center and the secondary radiation are calculated 
by radiation transfer throughout the slabs downstream. 

The dust-to-gas ratio (d/g) and the abundances of He, C, N, 
O, Ne, Mg, Si, S, A, Fe relative to H, are also accounted for. 
They afl'ect the calculation of the cooling rate. The dust grains 
are heated radiatively by photoionization and collisionally by 
the shock. The distribution of the grain radii downstream re- 
sults from sputtering. The geometrical thickness of the emitting 
nebula (D), determines whether the model is radiation-bound or 
matter-bound. 

The input parameters providing the best fit of the line ratios 
establish the physical conditions in each of the observed posi- 
tions. 

3.3. The code: calculation process 

The calculations initiate at the shock front where the gas is com- 
pressed and thermalized adiabatically, reaching a maximum tem- 
perature in the immediate post- shock region (T ~ 1.5 x 10^ 
(Vs /lOO kms"^)^). T decreases downstream following recom- 
bination. The cooling rate is calculated in each slab. The down- 
stream region is cut up into a maximum of 300 plane-parallel 
slabs with diff'erent geometrical widths calculated automatically, 
in order to account for the temperature gradient (Contini 1997 
and references therein). 

In each slab, compression is calculated by the Rankine- 
Hugoniot equations for the conservation of mass, momentum 
and energy throughout the shock front. Compression (n/no) 
downstream ranges between 4 (the adiabatic jump) and > 10, de- 
pending on Vs and ^o- The stronger the magnetic field, the lower 
the compression downstream is, while a higher shock velocity 
corresponds to a higher compression. 

The ionizing radiation from an external source is character- 
ized by its spectrum and by the flux intensity. The flux is cal- 
culated at 440 energies, from a few eV to keV. Due to radiative 
transfer, the spectrum changes throughout the downstream slabs, 
each of them contributing to the optical depth. In addition to the 
radiation from the primary source, the efl'ect of the difl'use radia- 
tion created by the gas line and continuum emission is also taken 
into account, using 240 energies to calculate the spectrum. 

For each slab of gas, the fractional abundance of the ions 
of each chemical element is obtained by solving the ioniza- 
tion equations. These equations account for the ionization mech- 
anisms (photoionization by the primary and difl'use radiation, 
and collisional ionization) and recombination mechanisms (ra- 
diative, dielectronic recombinations), as well as charge transfer 
efl'ects. The ionization equations are coupled to the energy equa- 
tion if collision processes dominate, and to the thermal balance 
if radiative processes dominate. The latter balances the heating 
of the gas due to the primary and diff'use radiations reaching the 



slab with the cooling due to recombinations and collisional ex- 
citation of the ions followed by line emission, dust collisional 
ionization and thermal bremsstrahlung. The coupled equations 
are solved for each slab, providing the physical conditions nec- 
essary for calculating the slab optical depth, as well as its line 
and continuum emissions. The slab contributions are integrated 
throughout the cloud. 

In particular, the absolute line fluxes corresponding to the 
ionization level / of element K are calculated by the term n]^(i), 
which represents the density of the ion /. We consider that 
^k(0 = X(0[K/H]^H, where X(/) is the fractional abundance 
of the ion i calculated by the ionization equations, [K/H] is the 
relative abundance of the element K to H and nu is the density 
of H (in number cm"^). In models including shock, nu is cal- 
culated by the compression equation (Cox 1972) in each slab 
downstream. So the abundances of the elements are given rela- 
tive to H as input parameters. 

Dust grains are coupled to the gas across the shock front by 
the magnetic field (Viegas & Contini 1994). They are heated by 
radiation from the AGN and collisionally by the gas to a max- 
imum temperature which is a function of the shock velocity, of 
the chemical composition and of the radius of the grains, up 
to the evaporation temperature (rdust ^ 1500 K). The grain ra- 
dius distribution downstream is determined by sputtering, which 
depends on the shock velocity and on the density. Throughout 
shock fronts and downstream, the grains might be destroyed by 
sputtering. 

Summarizing, the code starts by adopting an initial Tq (~ 10^ 
K) and the input parameters for the first slab (Sections 3.1 and 
3.2). It then calculates the density from the compression equa- 
tion, the fractional abundances of the ions from each level for 
each element, line emission, free-free emission and free-bound 
emission. It re-calculates Tq by thermal balancing or the enthalpy 
equation, and calculates the optical depth of the slab and the 
primary and secondary fluxes. Finally, it adopts the parameters 
found in slab / as initial conditions for slab /-hi. The line inten- 
sities and the bremsstrahlung are calculated from the gas which 
emits the spectrum. The data are observed at Earth, therefore 
they diverge by a factor which depends on the distance of the 
nebula from the active centre (r), and on the distance (d) of the 
galaxy to Earth (r^/d^). We then calculate the line ratios to a 
specific line (in the present case HjS which is a strong line), and 
we compare them with the observed line ratios. 

4. Modelling the NGC 7212 spectra 

The spectra are rich in number of lines, therefore the calculated 
spectra are strongly constrained by the observed line ratios. The 
line ratios are diff'erent in each of the observed spectra, revealing 
diff'erent physical conditions from region to region. At each stage 
of the modelling process, if a satisfactory fit is not found for all 
the lines, a new iteration is initiated with a diff'erent set of input 
parameters. When a ratio is not reproduced, we check how it de- 
pends on the physical parameters and decide consequently how 
to change them. Each ratio has a diff'erent weight, however we 
generally consider that the observed spectrum is satisfactorily 
fitted by a model when the strongest line ratios are reproduced 
by the calculation within a 20% discrepancy and weak line ratios 
within 50%. 

We have run a grid of models (~ 100 for each position) in 
order to select the most appropriate one. The fit was successful 
in 111 out of the 256 spectra, since they were the most com- 
plete in number of lines. Fig. [3] shows a map of the S/N ratio of 
the stellar continuum with the number of the modelled spectrum 
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Fig. 3. 2D map of the S/N ratio of the stellar continuum. The 
numbers indicate the spectra successfully modelled. 
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Fig. 4. The correlation of the input parameters with calculated 
[OIII]5007+/[OII]3727+ line ratios. 



overlaid. Most of these spectra have S/N larger than 6. The final 
gap between observed and calculated line ratios is due to both 
random and systematic observational errors, as well as to the 
uncertainties of the atomic parameters adopted by the code, such 
as recombination coefficients, collision strengths etc., which are 
continuously updated, and to the choice of the model. Models 
are generally allowed to reproduce the data within a factor of 2. 
This leads to input parameter ranges of a few per cent. 

4.1. First step: choice of physical parameters 

In the first iterations, we try to reproduce the 
[OIII]5007+4959/HyS line ratio (5007+4959 will be re- 
ferred to in the following as 5007+ ; the + indicates that the 
doublet 5007, 4959 is summed up), the highest ratio in general, 
by readjusting F and Vs (Vs however, is constrained throughout 



a small range by the observed FWHM). The higher F, the 
higher the [OIII]/HyS and the [OIII]/[OII] line ratios are (see 
Fig.©, as well as Hell/HyS. Moreover, a high F maintains the 
gas ionized far from the source, yielding enhanced [01] and 
[SII] lines. Recall that these lines behave similarly because 
S first ionization potential (10.36 eV) is lower than that of O 
(13.61 eV). 

Then, the [O II] 3727+ doublet is considered. If the flux from 
the active center is low (F < 10^ ph cm"^ s"^ eV"^), a shock 
dominated regime is found, which is characterised by relatively 
high [O II]/[0 III] ( > 1, as can be seen in Fig.©. These lines can 
be drastically reduced by collisional deexcitation at high electron 
densities (Hq > 3000 cm"^). 

The gas density is a crucial parameter. In each cloud, it 
reaches its upper limit downstream and remains nearly constant, 
while the electron density decreases following recombination. A 
high density, increasing the cooling rate, speeds up the recombi- 
nation process of the gas, enhancing the low ionization lines. In 
fact, each line is produced in a region of gas at a different Uq and 
Tq, depending on the ionization level and the atomic parameters 
characteristic of the ion. The density n, which can be roughly 
inferred from the [SII] 67 16/6731 doublet ratio, is related with 
no by compression downstream (n/no), which ranges between 
~ 4 and ~ 100, depending on Vs and ^o- The [SII] lines are 
also characterised by a relatively low critical density for colli- 
sional deexcitation. In NGC 7212 the [S II]6716/6731 line ratios 
are ~ 1 in all observed positions, constraining the choice of hq. 
In some cases the [SII] 67 16/6731 line ratio varies from >1 to 
< 1 throughout a relatively small region, since the [S II] line ra- 
tios depend on both the temperature and electron density of the 
emitting gas (Qsterbrock 1 1 989b . The median error of the mod- 
elled ratios is about 10%, and 90% of this sample has an error 
lower than 30%. 

We fit the lines of a single spectrum adopting an no which 
is combined with the other input parameters, including the pre- 
shock magnetic field ^o- The magnetic field has an important 
role in models accounting for the shock. The stronger Bq, the 
lower the compression is. So, lower densities are compensated 
by a lower magnetic field. If hq in one spectrum is inconsistent 
with no calculated in the surrounding regions, we can reduce the 
no divergence by changing the magnetic field, because no and 
^0 are strongly (anti)correlated. This procedure would lead to a 
smoother distribution of no throughout the ENLR of NGC 7212, 
which is however not expected because collisions and merg- 
ing phenomena render the density field dishomogeneous. So we 
have adopted a constant = 10"^ G, which is suitable for the 
ENLR of AGN (Beck 2011). This choice is slightly unreaHstic 
but it is justified by the fact that the preshock densities range 
between values within a factor of 2 (~ 80-150 cm"^, Table[T]). 
Table 1 shows that could be modified in regions 103-105, 
121-122, 132-133, and 156-157. These regions are too few as to 
yield a trend to throughout the ENLR. Moreover, ffuctua- 
tions would not affect the ensemble of the results. 

The final choice of a model is dictated by [OI]6300+/H/5 , 
which is calculated consistently with [O Ill]/Uj3 and [O ll]/Bfi. 
The geometrical thickness (D) of the emitting clouds is mainly 
determined by [0 1]/UJ3 because [O I] is emitted from gas at rel- 
atively low temperature, at the cloud edge opposite the shock 
front. 

A first analysis of the data shows that the models are mostly 
matter bound, namely, the integration of the calculated line 
fluxes downstream is interrupted at such a distance from the 
shock front that all the line ratios of a single spectrum are sat- 
isfactorily reproduced (by 10 per cent for the strongest ones). 
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Fig. 5. Top: histograms of the shock velocities (left), pre-shock density (middle) and the AGN ionizing flux (right). Bottom: 2D 
maps of the same parameters. The units are km s"^ (left), cm"^ (middle) and ph cm"^ s"^ eV"^ (right) 

. Spatial scale is 1 arcsec px~^ Contours of the stellar continuum are overlapped onto the map of the ionizing flux. North is up and East is to the 

left. 



The low [OI]6300/HyS (< 0.6, Fig.[2]) observed by Cracco et al. 
are hardly reproduced by radiation bound models, which are re- 
vealed in a few cases (represented by the maximum D in Fig.|4]). 
The [S II] 67 16, 673 1 lines are also emitted by relatively cool gas. 
Therefore, their ratios to HyS are determined by [O I]/HyS. 

So far, the results of modelling correspond to the set of input 
parameters presented in Table[T] The relative maps are shown 
in Fig. [5] The uncertainties of the input parameters are approxi- 
mately 10%. In fact, a change of 10% corresponds to variations 
of the resulting line ratios within 20% for the strongest lines up 
to 50% for the weakest lines. For all the models, solar abun- 
dances were used. In Table[2]we compare the spectra calculated 
in each position with the observed corrected line ratios to HyS. 
Each of the observed line ratios is followed by the calculation 
results. The numbers in the first column refer to Fig.[T] Not all 
the lines used for modelling are included in Table[2l We have fo- 
cused on [OIII]5007-h, [NeIII]3869-h, [OII]3727-h, [OI]6300-h, 
Hell 4686, Hel 5876, [NII]6548-h, [S II]6716,6731. For aU 
the spectra, we assumed that the reddening corrected ratio is 
HQf/HyS=2.8, model calculations give Ha/HyS ~ 2.9 and Hy/HyS 
~ 0.46. 

Tabled and Fig.[6] show that the observed [OIII]/HyS, 
[O II]/HyS and [O I]/Hy^5 are reproduced by the models with rel- 
atively high precision (the discrepancy is < 20 %). 

The Ne/H relative abundance used in the calculation of the 
spectra is 10""^. Comparison with the data shows that Ne/H 
higher by a factor of 1.5-2 should be adopted. This would ap- 
proac hthe N e/H = 1.95 x 10"^ solar value calculated by Bahcall 
et al. (120051) . Table[2]shows that [Fe VII]6087/HySis not negligi- 
ble in some regions, in the nucleus and in the regions to the East 



of the nucleus, even if the iron line is faint and the median error 
is about 0.3. 

The modelled He I/HyS and He II/HyS ratios are in agreement 
with the observed values, within their large error bars. [N II]/HyS , 
[S II]/H/5 line ratios show a large scatter (see Fig.O, which does 
not depend on the observational errors. We point out that at this 
stage of modelling, we considered solar abundances. We will 
now determine the metallicity relative to N and S, by constrain- 
ing in the following calculations the relative abundances of N/H 
and S/H, in order to obtain a perfect fit of calculated to observed 
line ratios. 



4.2. Second step: relative abundances 

The ratios of calculated to observed [SII]/HyS and [NII]/HyS in 
each position are shown in Tabled Unfortunately there are no 
data for S lines from higher ionization levels, which could indi- 
cate whether the choice of the model is misleading, or diff'erent 
relative abundances should be adopted. Since sulphur can be eas- 
ily depleted from the gaseous phase and trapped into dust grains 
and molecules, we considered that the [S II]/Hy5 line ratios indi- 
cate the S/H relative abundance. 

The same is valid for the [Nil] lines. The observations do 
not contain any other strong line of N, which could confirm 
whether the discrepancy between the calculated and the ob- 
served [N II]/HyS line ratios can be resolved by changing the N/H 
relative abundance or the other input physical parameters. Solar 
nitrogen is less abundant than other solar elements e.g. oxy- 
gen, neon, etc., therefore it is not an efficient coolant, namely, it 
does not aff'ect strongly the cooling rate of the gas downstream. 
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Table 1. Physical parameters resulting from modelling the observed spectra. 



pos 


77 — T 


T 

no 




pi 


— 7^ 

D"^ 


Uj3 abs^ 


pos 


— T 


T 

no 


775 

pi 


7771 


H/? abs^ 


35 


300. 


90. 


20 


7.8 


0.096 


123 


750. 


90 


16 


0.74 


0.04 


36 


300. 


100. 


70 


9.5 


0.25 


124 


60. 


60. 


1 


860 


0.0047 


37 


160. 


100. 


4 


0.76 


0.0073 


124 


600. 


80. 


30 


54 


0.3 


38 


110. 


110. 


9 


5.6 


0.02 


125 


600. 


90. 


30 


63 


0.39 


39 


120. 


110. 


8 


2.9 


0.0165 


132 


550. 


60. 


95 


50 


0.5 


40 


120. 


110. 


8 


3.5 


0.018 


133 


110. 


110. 


10 


6 


0.022 


51 


120. 


100. 


12 


230 


0.04 


134 


110. 


110. 


10 


6 


0.022 


52 


90. 


130. 


12 


17 


0.038 


135 


110. 


110. 


10 


6 


0.022 


53 


80. 


120. 


12 


15 


0.029 


136 


110. 


110. 


10 


13 


0.025 


54 


110. 


130. 


10 


5. 


0.032 


137 


100. 


110. 


10 


30 


0.026 


55 


140. 


130. 


9 


2.3 


0.026 


138 


100. 


120. 


12 


280 


0.04 


56 


110. 


110. 


10 


5.6 


0.022 


139 


100. 


120. 


10 


280 


0.042 


57 


110. 


100. 


10 


6.8 


0.02 


139 


750. 


90. 


13 


0.5 


0.057 


58 


180. 


100. 


18 


3.6 


0.025 


140 


90. 


90. 


11 


73 


0.046 


60 


750. 


110. 


80 


0.1 


0.13 


140 


750. 


80. 


13 


0.92 


0.037 


66 


160. 


100. 


15 


4. 


0.026 


141 


90. 


90. 


11 


1000 


0.06 


67 


130. 


100. 


5 


4.7 


0.015 


141 


600. 


90. 


20 


27.2 


0.156 


68 


110. 


100. 


4.5 


5. 


0.01 


142 


80. 


90. 


1 


2.4 


0.0023 


69 


110. 


100. 


4.5 


4.4 


0.01 


148 


380. 


120. 


33 


3.5 


0.1 


70 


110. 


110. 


10 


6. 


0.027 


149 


110. 


110. 


9 


7.6 


0.021 


71 


550. 


80. 


50 


9.5 


0.3 


150 


220. 


100. 


10 


1.6 


0.02 


72 


220. 


190. 


32 


1. 


0.09 


151 


160. 


100. 


5 


1.7 


0.011 


73 


220. 


100. 


8 


1.23 


0.02 


152 


160. 


100. 


4 


1.3 


0.009 


74 


220. 


110. 


10 


1.33 


0.025 


153 


160. 


100. 


4 


1.75 


0.01 


83 


60. 


100. 


11 


50 


0.019 


154 


70. 


90. 


0.7 


2 


0.0016 


84 


100. 


90. 


5 


170 


0.013 


155 


70. 


90. 


0.7 


1.6 


0.0015 


85 


110. 


100. 


11 


12 


0.03 


155 


750. 


90. 


1 


0.2 


0.0013 


85 


110. 


100. 


7 


5.5 


0.02 


156 


50. 


100. 


0.07 


0.15 


0.00009 


86 


320. 


100. 


60 


7 


0.23 


156 


600. 


60. 


0.04 


0.5 


0.0003 


87 


320. 


100. 


15 


7 


0.28 


157 


650. 


60. 


0.05 


0.5 


0.0003 


88 


280. 


120. 


15 


7 


0.28 


165 


110. 


110. 


13 


13 


0.039 


89 


120. 


150. 


13 


7 


0.04 


166 


220. 


100. 


8 


1.2 


0.019 


90 


80. 


90. 


0.3 


2.7 


0.0001 


167 


160. 


100. 


5 


1.4 


0.011 


91 


120. 


100. 


5 


5 


0.013 


168 


160. 


100. 


4 


1.2 


0.0087 


91 


800. 


100. 


13 


0.24 


0.028 


169 


160. 


100. 


3.5 


1.1 


0.0085 


92 


120. 


100. 


4 


4.7 


0.013 


170 


70. 


90. 


0.7 


1.7 


0.0016 


92 


700. 


100. 


13 


0.33 


0.039 


171 


70. 


100. 


0.25 


0.09 


0.000135 


93 


650. 


90. 


13 


0.49 


0.036 


172 


50. 


100. 


0.07 


0.12 


0.000067 


99 


120. 


110. 


8 


3.8 


0.018 


173 


80. 


70. 


0.05 


0.15 


0.0001 


100 


550. 


80. 


45 


10 


0.28 


181 


300. 


120. 


26 


1.7 


0.07 


101 


550. 


80. 


70 


22 


0.45 


182 


160. 


100. 


5 


1.4 


0.011 


102 


550. 


70. 


50 


8.3 


0.2 


183 


120. 


100. 


3 


1.5 


0.006 


103 


550. 


70. 


50 


12.4 


0.22 


184 


160. 


100. 


4 


1.2 


0.001 


104 


380. 


120. 


33 


1.5 


0.11 


185 


80. 


90. 


1 


1.6 


0.0022 


105 


180. 


180. 


20 


1. 


0.05 


186 


70. 


90. 


0.7 


1.5 


0.0015 


106 


110. 


110. 


10 


13 


0.04 


187 


180. 


70. 


0.03 


0.2 


0.00019 


107 


50. 


100. 


0.5 


2.2 


0.001 


188 


180. 


70. 


0.06 


0.14 


0.00011 


107 


800. 


100. 


3 


0.156 


0.002 


197 


110. 


110. 


9 


6.1 


0.021 


108 


100. 


110. 


9 


280 


0.033 


198 


180. 


160. 


20 


1.6 


0.051 


108 


750. 


90. 


16 


0.58 


0.07 


199 


90. 


100. 


1.6 


1.2 


0.003 


109 


750. 


90. 


16 


0.58 


0.07 


200 


70. 


90. 


0.9 


2.7 


0.0018 


110 


750. 


90. 


13 


0.65 


0.06 


201 


70. 


90. 


0.7 


1.3 


0.0015 


111 


750. 


90. 


10 


0.4 


0.04 


202 


100. 


100. 


0.8 


0.4 


0.0017 


115 


600. 


80. 


65 


12 


0.36 


203 


160. 


60. 


0.9 


1 


0.002 


116 


550. 


80. 


60 


7 


0.28 


213 


300. 


120. 


30 


1.8 


0.05 


117 


550. 


80. 


50 


9.5 


0.3 


214 


500. 


80. 


40 


7 


0.225 


118 


550. 


80. 


60 


11.3 


0.3 


215 


120. 


150. 


13 


2.7 


0.037 


119 


550. 


80. 


70 


9.2 


0.32 


216 


160. 


100. 


3.5 


1 


0.0085 


120 


380. 


120. 


33 


1.5 


0.11 


217 


90. 


90. 


11 


45 


0.0045 


121 


120. 


150. 


13 


4 


0.04 


218 


100. 


80. 


7 


20 


0.023 


122 


130. 


80. 


6 


37 


0.028 


219 


90. 


120. 


0.75 


0.73 


0.002 


123 


110. 


110. 


10 


15 


0.04 















Notes. in kms ^; in cm ^; in 10^ photons cm ^ s ^ eV ^ at the Lyman limit; '^'^^ in 10^^ cm; in erg cm ^ s ^ 
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Fig. 6. Comparison between the observed and the calculated emission line ratios for the lines in Table 2. cr is the median error for 
the observed ratios. The solid line is the 1:1 line. 



Varying N/H will change only the [N II]/HyS line ratio, but not the 
other line ratios, while changing 0/H the whole input parameter 
set must be readjusted. 



4.2.1. [Nil] and [ON] 

Tableland Fig.Oshow that the lines which systematically dis- 
agree with the observations are the [Nil] doublet at 6548+6584 
A. The doublet is summed up in the table. In the calculations, 
solar N/H value 9.1 x 10"^ (Allen 1976) was adopted. Table[2l 
shows that the observed [NII]/HyS line ratios are mostly over- 
predicted by the models. The discrepancies between observed 
and calculated values are shown in Fig.[71 [Nil] lines come from 



the same region of the gas that emits the [O II] lines. [O II]/HyS 
line ratios depend on the physical parameters and are constrained 
by [OIII]/HyS and [OI]/HyS, so the discrepancy in the [NII]/HyS 
line ratios can be explained by abnormal N/H relative abun- 
dances. In the left diagram of Fig.[7]([NII]/H/5)caic/([NII]/HyS)obs 
and ([O II]/HyS)caic/([0 II]/HyS)obs are plotted in the different posi- 
tions. If discrepancies between calculation and observation were 
due to the model, both [NII]/H/5 and [O II]/H/5 would show the 
same variations. Contrariwise, in positions beyond 100 (regions 
to the East of the nucleus) the [NII]/HyS calculated line ratios 
overpredict the data by a factor of at least 1.5, whereas the ob- 
served [O II]/HyS are always well reproduced by the models. This 
indicates that N/H should be depleted by a factor of at least 1.5, 
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# * 

* % * *# 

* # # ***** 


* % * * ** 
* ** ^ f ** ^ * 


* * 






Fig. 7. Variations of calculated/observed line ratios as a function of the number of spectra (see Fig.[T]). Left: [NII]/HyS (red aster- 
isks), [OII]/HyS (blue asterisks); middle : [OI]/HyS (green asterisks), [S II]67 16+673 1/HyS (black asterisks); right: Hell/HyS (magenta 
asterisks). 






0.48 0.63 0.66 0.7 0.76 0.81 0.86 0.93 0.99 1.1 1.2 1.4 2. J 0.42 0.47 0.52 0.57 0.62 0.66 0.71 0.76 0.81 0.86 0.9 0.95 1 ) . 27 1 . 7 5 . 8 82 . 9 2 7 . 9 6 5 1 . 2 1.1 1.17 1.24 1.35 1.56 1.78 3.4 

Fig. 8. The obs/calc maps of nitrogen (left), oxygen (middle) and sulphur (right). Values lower than 1 indicate depletion. Spatial 
scale is 1 arcsec px"^. North is up and East is to the left. 



while 0/H is almost solar. The map showing the N/H depletion 
and excess (TableO appears in Fig.[8l 



4.2.2. [Sll] and [01] 

From the last two columns of Table[2] we can realize that the dou- 
blet ratios [S II] 67 16/673 1 ~ 1 , indicative of the emitting gas den- 
sity, are relatively low (< 10^ cm"^). The [S II] and [O I] lines are 
emitted from relatively cool gas in similar physical conditions. 
In Fig.[7] (middle diagram) both ([SII]/HyS)caic/([SII]/HyS)obs and 
([OI]/HyS)caic/([OI]/HyS)obs ratios show up. Table[51and Fig.[6l 
suggest that [SII] line ratios to HjS calculated adopting solar 
abundances are acceptable, except in some positions. At the 
North-East of the nucleus the calculated [S II]/HyS ratios overpre- 
dict the data and at other positions they even underpredict them 
(see Fig. [7]). The map of the S/H depletion and excess (TableO 
is shown in Fig.[8l Eventually, the [S II] observed lines include 
a large contribution from the ISM cool gas. The [SII]4070+ 
lines are not included in Tabled because they do not appear 
in all the spectra and since they are blended with other lines 
at nearby wavelengths. However, model results indicate that the 
[S II]4070+/H/5 ratios are relatively high for high velocities and 
reproduce the data in a few positions. 

4.2.3. He II and the oxygen lines 

The first cycle of calculations yields He ll/H/S line ratios gener- 
ally lower than those observed. Indeed, higher He 11/UJ3 can be 



obtained by increasing the photoionization flux, which however 
spoils the goodness of fit of the [OIII]/HyS and [OII]/HyS line 
ratios. Another option is to reduce the geometrical thickness of 
some of the clouds, excluding the cooling region. Accounting 
only for the high temperature zone close to the shock front 
downstream, only the Hell lines would be strong. Summing up 
the spectra of these thinner clouds with the other ones would 
give high Hell/HyS ratios. The three most important line ratios 
[O III]/HyS, [O II]/HyS and [O I]/HyS are all well fitted adopting the 
solar 0/H relative abundance, 6.6 x 10""^ (Allen I 19761) . However, 
Cracco et al. (2011), compared the line ratios with models cal- 
culated by pure photoionization and found a low metallicity in 
NGC 7212. Oxygen is a strong coolant. Decreasing 0/H in the 
models will induce a lower cooling rate downstream, yielding 
a larger zone of hot gas capable of emitting higher He II/UJS. 
Moreover, the Hell line sits on top of the [Felll] multiplet, 
whose most prominent lines are at 4658 A and 4702 A. The lines 
are blended with Hell in high velocity regimes. Consequently, 
we adopted in the next cycle of calculations a lower 0/H rela- 
tive abundance in positions where the observed He 11/UJ3 lines 
ratios were underpredicted. This required an adjustment of the 
input parameters also regarding the physical conditions. Finally, 
we added to He 11/UJ3 the [Fe III]/HyS line ratios, which can be 
as high as He 11/UJ3 in some positions. The comparison of calcu- 
lated with observed Hell/HyS improved (Fig.[7l right diagram) 
except in some positions, e.g. 109, 124, 200, where observa- 
tional errors are high (larger than 60%). In particular, we ob- 
tained 0/H = 3.6 X 10-4 in positions 60 and 115, 4.6 x lO'^ in 
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16 16.5 17 17.5 18 18.5 19 19.5 

log (distance from shock front) 



Fig. 9. The profile of HI (solid line), Hel (dot-dashed) and Hell 
(dashed) fractional abundances as a function of distance from 
the shock front for position 53. 



positions 116, 141, 187, 188, 203 and 4 x 10"^ in 213. Moreover, 
0/H = 3 X 10-4 was used in positions 125, 140, and 2.6 x 10""^ 
in position 133. Interestingly, most of these positions refer to 
relatively high shock velocities (see Table[T]). For all other posi- 
tions, the solar 0/H value was adopted. The map of the oxygen 
depletion is shown in Fig.[8l 

4.2.4. Hel/HyS line ratios 

The Hel/HyS diagram in Fig. [6] and the results in Table|2] show 
that the calculated Hel/HyS line ratios mostly range between 0.11 
and 0.13, while the observations range between 0.04 and 0.4. 
To understand the results, we show in Fig. [9] the profile of the 
temperature and of HI, Hel and Hell fractional abundances as a 
function of distance from the shock front, which corresponds to 
the cloud edge reached by the flux from the AGN. We have cho- 
sen position 53 which shows a strong F and a relatively low Vs . 
The temperature maximum between 10^^ and 3x10^^ cm is real 
and not due to the calculation approximation. It arises when hy- 
drogen becomes neutral, followed by the decrease of the density 
of electrons, which, colliding with the heavy elements, produce 
cooling downstream (Williams 1967). The calculations of the He 
ion fractional abundances follow the temperature fluctuations by 
small variations, while the HI fluctuation is imperceptible, so the 
Hel/HyS calculated range is small. On the other hand, the obser- 
vational error of Hel reaches 79%. 

4.2.5. [Fell] 

Rodriguez- Ardila et al. (2004) studied the [Fe II] line emission 
in AGN and found that [FeII]/Pa/5 increases from pure pho- 
toionization to a shock excitation regime. For instance, [Fell] 
emission and double profiles are due to jet-gas interaction in 
Mrk 78 (Whittle et al. [T9881 Ramos Almeida et al. 2006) and 
[Fell] lines are strong in the radio-jet interaction region in Mrk 
34 (Jackson & Beswick ,2007i) . This could indicate the impor- 
tance of the shocks for strong [Fell] lines. Iron first ionization 
potential (7.9 eV) which corresponds to the [Fell] lines is lower 
than that of O^ and H^. Therefore [Fell] lines are mainly emit- 
ted from cool recombined gas. Collisional ionization coefficients 
are low at those low temperatures (< 10"^ K), so strong [Fell] 
lines in the shock-dominated case (as for supernovae) are due to 



avy elements in NGC 7212 

Table 3. Calculated relative abundances. 



pos. 


N/H 




S/H 




pos. 


N/H 




S/H 




35. 


1.17x10- 


-4 


6.99x10- 


-6 


123. 


1.08x10- 


-4 


1.32x10" 


-5 


36. 


1.08x10- 


-4 


1.27x10- 


-5 


124. 


1.01x10- 


-4 


1.67x10" 


-5 


37. 


6.20x10- 


-5 


1.20x10- 


-5 


124. 


1.56x10- 


-4 


2.09x10" 


-5 


38. 


7.19x10- 


-5 


1.27x10- 


-5 


125. 


2.28x10- 


-4 


1.47x10" 


-5 


39. 


7.22x10- 


-5 


1.51x10- 


-5 


132. 


1.05x10- 


-4 


1.48x10" 


-5 


40. 


6.81x10- 


-5 


1.34x10- 


-5 


133. 


8.64x10- 


-5 


1.07x10" 


-5 


51. 


1.05x10- 


-4 


1.08x10- 


-5 


134. 


8.62x10- 


-5 


1.07x10" 


-5 


52. 


1.21x10- 


-4 


1.36x10- 


-5 


135. 


8.75x10- 


-5 


1.01x10" 


-5 


53. 


6.96x10- 


-5 


1.55x10- 


-5 


136. 


9.02x10- 


-5 


8.95x10" 


-6 


54. 


8.08x10- 


-5 


1.69x10- 


-5 


137. 


8.94x10- 


-5 


8.31x10" 


-6 


55. 


5.85x10- 


-5 


1.70x10- 


-5 


138. 


1.14x10- 


-4 


8.34x10" 


-6 


56. 


5.96x10- 


-5 


1.84x10- 


-5 


139. 


6.62x10- 


-5 


1.70x10" 


-5 


57. 


6.08x10- 


-5 


1.69x10- 


-5 


139. 


9.75x10- 


-5 


1.19x10" 


-5 


58. 


4.36x10- 


-5 


1.47x10- 


-5 


140. 


3.05x10- 


-5 


1.33x10" 


-5 


60. 


1.09x10- 


-4 


2.77x10- 


-5 


140. 


5.86x10- 


-5 


1.44x10" 


-5 


66. 


7.79x10- 


-5 


1.77x10- 


-5 


141. 


1.06x10- 


-4 


1.17x10" 


-5 


67. 


1.10x10- 


-4 


1.26x10- 


-5 


141. 


1.06x10- 


-4 


7.51x10" 


-6 


68. 


7.88x10- 


-5 


1.25x10- 


-5 


142. 


6.77x10- 


-5 


6.56x10" 


-6 


69. 


6.80x10- 


-5 


1.31x10- 


-5 


148. 


8.53x10- 


-5 


1.59x10" 


-5 


70. 


7.71x10- 


-5 


2.21x10- 


-5 


149. 


8.35x10- 


-5 


1.17x10" 


-5 


71. 


1.41x10- 


-4 


2.54x10- 


-5 


150. 


7.62x10- 


-5 


1.07x10" 


-5 


72. 


5.52x10- 


-5 


2.00x10- 


-5 


151. 


7.40x10- 


-5 


8.50x10" 


-6 


73. 


5.46x10- 


-5 


1.69x10- 


-5 


152. 


7.28x10- 


-5 


8.07x10" 


-6 


74. 


5.98x10- 


-5 


1.45x10- 


-5 


153. 


7.54x10- 


-5 


7.17x10" 


-6 


83. 


6.12x10- 


-5 


1.67x10- 


-5 


154. 


6.91x10- 


-5 


5.70x10" 


-6 


84. 


8.82x10- 


-5 


7.37x10- 


-6 


155. 


6.07x10- 


-5 


7.68x10" 


-6 


85. 


7.84x10- 


-5 


1.29x10- 


-5 


155. 


5.93x10- 


-5 


7.79x10" 


-6 


86. 


8.85x10- 


-5 


2.14x10- 


-5 


156. 


5.90x10- 


-5 


6.51x10" 


-6 


87. 


8.40x10- 


-5 


2.33x10- 


-5 


156. 


5.46x10- 


-5 


7.99x10" 


-6 


88. 


8.05x10- 


-5 


2.51x10- 


-5 


157. 


5.66x10- 


-5 


8.12x10" 


-6 


89. 


7.00x10- 


-5 


1.74x10- 


-5 


165. 


9.13x10- 


-5 


1.27x10" 


-5 


90. 


8.47x10- 


-5 


1.30x10- 


-5 


166. 


7.33x10- 


-5 


1.07x10" 


-5 


91. 


6.42x10- 


-5 


2.80x10- 


-5 


167. 


6.56x10- 


-5 


1.01x10" 


-5 


91. 


7.46x10- 


-5 


1.19x10- 


-5 


168. 


6.44x10- 


-5 


9.05x10" 


-6 


92. 


5.38x10- 


-5 


1.19x10- 


-5 


169. 


6.35x10- 


-5 


8.05x10" 


-6 


92. 


8.11x10- 


-5 


9.13x10- 


-6 


170. 


6.35x10- 


-5 


6.37x10" 


-6 


93. 


5.04x10- 


-5 


1.30x10- 


-5 


171. 


7.80x10- 


-5 


8.35x10" 


-6 


99. 


9.07x10- 


-5 


1.42x10- 


-5 


172. 


5.95x10- 


-5 


6.56x10" 


-6 


100. 


1.52x10- 


-4 


1.89x10- 


-5 


173. 


1.05x10- 


-4 


5.68x10" 


-6 


101. 


1.54x10- 


-4 


1.74x10- 


-5 


181. 


9.43x10- 


-5 


1.79x10" 


-5 


102. 


1.11x10- 


-4 


1.73x10- 


-5 


182. 


6.80x10- 


-5 


1.12x10" 


-5 


103. 


1.16x10- 


-4 


1.93x10- 


-5 


183. 


5.85x10- 


-5 


1.11x10" 


-5 


104. 


8.60x10- 


-5 


1.73x10- 


-5 


184. 


6.13x10- 


-5 


1.02x10" 


-5 


105. 


7.28x10- 


-5 


1.43x10- 


-5 


185. 


5.55x10- 


-5 


8.35x10" 


-6 


106. 


1.08x10- 


-4 


1.43x10- 


-5 


186. 


5.72x10- 


-5 


7.67x10" 


-6 


107. 


7.23x10- 


-5 


1.68x10- 


-5 


187. 


6.81x10- 


-5 


3.74x10" 


-6 


107. 


6.72x10- 


-5 


9.40x10- 


-6 


188. 


6.96x10- 


-5 


4.33x10" 


-6 


108. 


6.71x10- 


-5 


1.20x10- 


-5 


197. 


7.82x10- 


-5 


1.66x10" 


-5 


108. 


9.21x10- 


-5 


9.25x10- 


-6 


198. 


8.02x10- 


-5 


1.44x10" 


-5 


109. 


9.89x10- 


-5 


1.36x10- 


-5 


199. 


5.36x10- 


-5 


1.19x10" 


-5 


110. 


8.17x10- 


-5 


9.04x10- 


-6 


200. 


5.21x10- 


-5 


8.33x10" 


-6 


111. 


7.00x10- 


-5 


1.38x10- 


-5 


201. 


5.96x10- 


-5 


1.04x10" 


-5 


115. 


9.10x10- 


-5 


1.43x10- 


-5 


202. 


6.01x10- 


-5 


9.77x10" 


-6 


116. 


1.31x10- 


-4 


1.43x10- 


-5 


203. 


6.13x10- 


-5 


5.78x10" 


-6 


117. 


1.46x10- 


-4 


1.57x10- 


-5 


213. 


6.37x10- 


-5 


2.03x10" 


-5 


118. 


1.36x10- 


-4 


1.77x10- 


-5 


214. 


1.43x10- 


-4 


2.14x10" 


-5 


119. 


1.29x10- 


-4 


1.79x10- 


-5 


215. 


8.90x10- 


-5 


1.53x10" 


-5 


120. 


9.38x10- 


-5 


1.52x10- 


-5 


216. 


5.76x10- 


-5 


9.79x10" 


-6 


121. 


1.00x10- 


-4 


1.07x10- 


-5 


217. 


9.17x10- 


-5 


1.39x10" 


-5 


122. 


1.45x10- 


-4 


1.30x10- 


-5 


218. 


9.57x10- 


-5 


1.53x10" 


-5 


123. 


7.28x10- 


-5 


1.21x10- 


-5 


219. 


6.61x10- 


-5 


6.97x10" 


-6 



Notes. N/H(solar)=9.1xlO-5;S/H(solar)= 1.6x10"^ 
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diffuse secondary radiation, which in turn is high as it is emitted 
from shock-heated gas slabs. 

In HIT regions, on the other hand, the black body flux is 
rapidly reduced by radiation transfer throughout the gas and the 
temperature drop, so the [Fell] lines are weak. In Fig.fTOl (top 
diagram) we illustrate the temperature and electron density pro- 
files downstream in positions 140 (left) and 149 (right). We have 
chosen such positions since they have very diff'erent shock veloc- 
ities and are reached by diff'erent flux intensities. Nevertheless, 
the spectra are sufl&ciently well reproduced by the models. The 
0^/0, 0^/0 and O^^/O fractional abundances and the frac- 
tional abundances of significant Fe ions are described, respec- 
tively, in the middle and the bottom diagrams. 

5. The continuum SED 

We checked the models by comparing them with the observed 
spectral energy distribution (SED) of the continuum, shown in 
Fig.[TT] The data (represented by black crosses) were extracted 
from the NASA/IPAC extragalactic database (NED) and are de- 
tailed in Table 4. They are taken with diff'erent apertures, and 
therefore they cover diff'erent regions of the galaxy. However, for 
apertures of few arcsec, the emission is very likely dominated by 
the AGN. We expected the same, in case of the hard X-ray emis- 
sion, even if the aperture is large. Some problems can arise in the 
soft X-ray and infrared domains, where the fluxes are emitted by 
a large fraction of the galaxy and additional contributions, e.g. 
supernova remnants and dust, could not be negligible. 

Most of the AGN emission in the optical-UV range is hidden 
by the molecular dusty torus, which makes visible the black- 
body flux from the old stellar population, here reproduced by 
black-body radiation at T = 5000 K. The archival GALEX data 
in the UV follow the black-body old star radiation. Nevertheless, 
the power-law flux from the AGN and the shocks heat the gas 
which contributes to the SED by its reprocessed bremsstrahlung 
emission. 

Fig.[TT] shows that the data cover the X-ray domain, 
indicating the presence of high temperature regions. The 
bremsstrahlung at high frequencies depends on the temperature 
of the gas in the downstream region behind the shock front. The 
higher the shock velocity, the higher the temperature of the gas 
downstream. Thus, the bremsstrahlung peaks at high frequencies 
for models calculated with high . To reproduce soft X-ray 
emission (0.5-2 keV) we have adopted the result of the model 
calculated with a shock velocity Vs = 550 km s"\ corresponding 
to position 71. NGC 7212 is a Compton thick source (Bennert 
et al. 2006) with reflected hard X-rays (2-10 keV) coming di- 
rectly from the AGN. Yet, hard X-rays could be contributed by 
shocks corresponding to Vs = 800 kms"^ as found in position 
107. Models with lower velocities do not contribute to the con- 
tinuum SED at high frequencies. 

The IR shows a composite black-body radiation from dust 
grains at difl'erent temperatures. Dust grains are heated radia- 
tively by the flux from the active center and collisionally by the 
gas. The data in the near IR domain are reproduced by dust re- 
processed radiation downstream of the highest velocity models 
adopting a dust-to-gas ratio d/g ~ 2 x 10""^ and 4 x 10"^ (by 
mass) in positions 71 and 107, respectively. We stress, for com- 
parison, that in the ISM d/g - 4x 10"^-4 x 10""^ (Contini & 
Contini 2007). 

The data in the radio range show the characteristic power-law 
of the synchrotron radiation created by the Fermi mechanism at 
the shock front. The free-free and free-bound continuum radia- 
tion flux are calculated consistently with the line emission fluxes 



-9 




log V (Hz) 

Fig. 11. Comparison of the continuum SED with the model re- 
sults. Black crosses: data from NED; black encircled crosses : 
the data from the NED believed to be dominated by the AGN; 
red asterisks: data from Sargsyan (2011); black squares: data 
from GALEX archive; green: synchrotron radiation; red: the 
dust reprocessed radiation flux calculated by models 71 (solid 
line) and 107 (dashed line); blue: the bremsstrahlung for position 
71 (solid line) and for position 107 (dashed line). Black line: the 
black body flux from the background stars. 

for each model. The bremsstrahlung emitted from the gas has the 
same frequency distribution in the IR-radio range for difl'erent 
models (see Contini et al. 2004) when free-free self-absorption 
is negligible (at rather low densities). 

6. Discussion 

6.1. Comparison with other observations 

To confirm the ranges of the parameters obtained by reproducing 
Cracco et al. (2011) data, we present in Table[5]the results from 
the modelling of Bennert et al. (2006) nuclear optical spectrum, 
corresponding to an aperture of V. We find a shock velocity 
of 370 km s"^ , slightly higher than that inferred from the aver- 
age observed FWHM (250 kms"^) provided by these authors. 
Near-IR spectroscopic data are reported by Ramos Almeida et 
al. ( 2009) for five Seyfert galaxies showing jet-gas interaction, 
including NGC 7212. Its IR nuclear spectrum (1.5 ''-aperture) is 
presented in Table[6l The line ratios are well reproduced with 
a model very similar to those which reproduce Cracco et al. 
( 201 1 ) optical spectra. The contemporary fit of neutral and coro- 
nal lines in diff'erent ranges requires composite models (see, e.g. 
Rodriguez- Ardfla et al. 2005 ). 

The shock velocities adopted to fit the Bennert et al. and 
Ramos Almeida et al. nuclear spectra, were compared to the val- 
ues obtained by the present modelling of the Cracco et al. spectra 
in a region of 3'' x 3'' centred on the nucleus. This region was 
chosen to take into account possible eff'ects of slit misalignen- 
ment with respect to the active nucleus. The values are similar: 
370, 500, and 110-550 kms"\ respectively. The densities are: 
190, 200, and 70-120 cm"^ respectively. RecaU that the densi- 
ties in the present modelling are preshock, corresponding to den- 
sities downstream higher by a factor of 4-10, depending on the 
shock velocity. The fluxes are also similar: logF = 10.85, 10, 
and 10-10.8, respectively. A difl'erence of factor -100 is appar- 
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log distance from shockfront [cm] log distance from shockfront [cm] 



Fig. 10. pos. 140 (left), pos. 149 (right); top diagram: the electron temperature (solid line) and electron density (dashed line) down- 
stream; middle diagram: the fractional abundances of (green), (black), O^^ (red); bottom diagram: the fractional abundances 
of Fe^ (black), Fe^^ (red), Fe^^ (magenta), Fe"^^ (cyan), Fe^^ (yellow), Fe^^ (blue). 



ent in the geometrical thickness of the clouds. Those referring to 
IR spectra by Ramos Almeida et al. are approximately 3x10^^ 
cm, while the other are about 2 x 10^^ cm. It seems that observ- 
able coronal lines are emitted from very thin clouds, in order to 
prevent recombination of the Si, S, etc. 

6.2. The physical picture of the NLR 

The 2D maps in Fig. [5] depict the distribution of the shock ve- 
locity Vs (left panel), of the preshock density no (central panel) 
and of the flux from the active nucleus F (right panel). They re- 
fer to the shock regime. High velocity shocks (600-800 km s"^) 
delineate the region where matter external to the ENLR collides 
with the ENLR gas. The highest velocities are found in the east- 
ern regions, where the maps show low no values. Matter is faster 
throughout the less dense gas regions, as expected from the mass 
conservation equation at the shock front. This confirms the pres- 
ence of shocks. The preshock densities demonstrate a narrow 
maximum near the eastern edge, which most probably obstruct 
the spreading of the high velocity shock. The histograms in 
Fig. [5] illustrate the distributions of the values, Vg , no and log F. 

Explaining the ionization of the perturbed gas in radio galax- 
ies, Tadhunter et al. (2000 ) discuss why the perturbed gas varies 
significantly in surface luminosity from object to object, while 
the quiescent gas does not. The correlation between F and Hfi 
absolute flux calculated at the nebula is shown in Fig. [121 as well 
as the correlation of the shock velocity with H/S. Hfi is chosen 
so as to represent the surface luminosity at the emitting nebula. 
Extrapolating from many difl'erent objects to difl'erent regions 
in the same galaxy. Fig. [12] shows that the luminosity of each 
region in NGC 7212 is strongly linked with the AGN flux reach- 
ing it. The H/S line emission is strong from regions at tempera- 
tures T >10'^ K. The correlation with is less straightforward. 
Indeed, at temperatures higher than 10^ K collisional parame- 
ters are much more important than the radiation ones, however, 
a higher generates higher densities downstream and conse- 
quently a higher cooling rate, which results in a reduced H re- 
combination zone. So we can hardly recognize the shocks result- 
ing from collisions by means of the object luminosity, but they 
can be found out by analysing line ratios and profiles. 
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Fig. 12. Shock velocities vs. HjS calculated flux (top panel); 
fluxes from the active centre vs. UJ3 calculated flux (bottom 
panel). Red: models calculated with high (> 700 km s"^); 
blue: the other models. 



The shock velocities calculated in the observed positions 
throughout the NGC 7212 ENLR, roughly show the regions 
where collision of matter from the ENLR with matter from out- 
side the cone occurred. We have found that the velocities are 
lower than 200 km s"\ and the preshock densities are relatively 
low (~ 100 cm"^) in all the positions representing the ENLR 
matter uncontaminated by high density debris from the BLR. 
These velocities refer to a smooth regime. Tadhunter et al. (120001) 
use the term "quiescent" to mean "in a quasi-stable dynamical 
configuration in a galaxy", that is, within a factor of 2 of the ro- 
tational velocity. The high velocities, on the other hand, refer to 
shocks. 

We confirm that shocks are strong beyond the edges of the 
cone overlying in Fig. [13] the shock velocities upon the contour 
map of the photoionization cone. Our calculations reveal that 
high velocity shocks are located at the edge of the ionization 
cone and beyond it, in agreement with the large FWHM of the 
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Table 4. List of photometric data extracted from NED. 



band 


logy 


log y Fv 


Instrument 


Aperture 


Kei. 


J— iU KeV 


io.Zj / 


-iZ.Zji 


AiVliVl 


jU 


i 


z— iu KeV 


1 Q 1 ^1 
io. iOi 


111^1 
-iZ. ioi 


AMM 


jU 


i 


Z— iU KeV 


1 O 1^1 

io.ioi 


1 o /ino 
-iZ.4Uo 


Chandra 


1 

i.J 


Z 


Z— iU KeV 


1 O 1/^1 

io.ioi 


-iZ.Zo/ 


Ciianara 


J.J X Z 


Z 


Z— iU KeV 


1 O 1^1 

io.ioi 


1 O 1 c? 
-iZ.io / 


Chandra 


i 


Z 


Z— iU KeV 


1 O 1/^1 

io.ioi 


-iZ.UoO 




i 


a 
J 


U. /— / KeV 


i /.yoo 


1 o o< 1 
-iZ.Zji 


A cr^ A 


i 


J 


u. /— z KeV 


i /.Ji J 


-i J.U>'D 


AoCA 


i 


J 


U.J— z KeV 


i /.4oU 


-i J.U4J 


VA /TA/r 

AiVliVl 


ju 


1 
i 


U.J— z Key 


1 T /I QH 

i /.4oU 


-i J.i /j 


Chandra 


J.J X Z 


Z 


r U V 


i j.zy 


-iz.zyy 


LrAi^HA 


/I T' 

4.Z 




JN U V 


1 < 1 1 o 

i j.i iZ 


1 1 TCQfl 
-i i. /JJ 




J.J 




r jjUW 


i4.yjo 


1 1 A.C\Aa 
-i i.D>'4 


rib 1/ACo 


U.J 


/I 
4 


r jjUW 


1 /I n<o 
i4.yjo 


1 1 '2<'7 
-i i.JJ / 


T_rcnr/ a /^c 
rib 1 /ACo 


1 " 
i 


/I 
4 


F330W 


14.958 


-11.054 




2.54" 


4 


nip 


14.833 


-10.591 




total mag 


5 




14.139 


-10.71 


2MASS 


Bulge mag 


6 




14.139 


-11.024^ 


2MASS 


AGN mag 


6 


lljjim 


13.397 


-10.3ir 


IRAS 


30" 


7 


25 jim 


13.079 


-10.035^ 


IRAS 


r 


7 


60 fim 


12.698 


-9.841^ 


IRAS 


r 


7 


100 yum 


12.477 


-9.832^ 


IRAS 


T 


7 


4.85 GHz 


9.685 


-14.623^ 


NRAO 


3.5' 


8 


4.85 GHz 


9.685 


-14.651 


NRAO 


3.7' X 3.3' 


9 


1.40 GHz 


9.146 


-14.747^ 


NRAO 


700" 


10 


365 MHz 


8.562 


-14.869^ 


UTRAO 


30" 


11 


73.8 MHz 


7.868 


-15.361^ 


VLA 


80" 


12 



References. (1) Guainazzi et al. (2005); (2) Levenson et al. (2006); (3) 
Ueda et al. (2001); (4) Munoz Marm et al. (2007); (5) Zwicky et al. 
(1965); (6) Peng et al. (2006); (7) Moshir et al. (1990); (8) Becker et 
al. (1991); (9) Gregory & Condon (1991); (10) White & Becker (1992); 
(11) Douglas et al. (1996); (12) Cohen et al. (2007). 



Notes. data selected to be representative of the AGN emission. 



Table 5. Comparison between our model and the Bennert et al. 
(2006) nuclear optical spectrum. 



line ratios to H/? 


obs. 


calc. 


[Oil] 3727+ 


2.63 


3.1 


[Nelll] 3869+ 


1.7 


1.3 


[OIII] 4363 


0.23 


0.1 


Hell 4686 


0.24 


0.2 


[OIII] 5007+ 


12.06 


12.02 


[FeVII] 6087 


0.04 


0.02 


[01] 6300 


0.56 


0.7 


[FeX] 6375 


0.03 


0.001 




2.87 


2.9 


[Nil] 6583 


1.99 


2.9 


[SII] 6716 


0.56 


0.3 


[SII] 6731 


0.63 


0.45 


^0 (cm"^) 




190 


Vs(kms-i) 




370 


5o ( gauss) 




10-4 


logF^ 




10.85 


D(cm) 




2 X 10^^ 


H/? (ergcm~^s~^) 




0.23 


n (cm~^) 




6330 



Notes. F is in photons cm ^ s ^ eV ^ at the Lyman limit. 



Table 6. Comparison between our model and the Ramos 
Almeida et al. (.2009J IR spectrum. 



line ratios to [S iii] 0.9 


obs. 


calc. 


[CI] 0.98 


0.02 


0.02 


[SVIII]0.99 


0.012 


0.01 


[SII] 1.03 


0.054 


0.045 


He 1 1.08 


0.35 


0.36 


[SIX] 1.25 


0.0057 


0.06 


[SIX] 1.43 


0.011 


0.012 


P« 1.96 


0.158 


0.2 


[SXI] 1.9 


0.003 


0.01 


[Si VI] 1.96 


0.028 


0.045 


riQ (cm"^) 




200 


Vs(kms-i) 




500 


Bq (gauss) 




10-4 


logF^ 




10. 


D (cm) 




3 X 10^5 


H/?(ergcm-2 s'O 




6.3 X 10-3 


H/?/[SIII] 




1.96 


n (cm~^) 




1300 



Notes. F is in photons cm ^ s ^ eV ^ at the Lyman limit. 




75 100 114 137 147 162 192 305 511 582 731 



Fig. 13. 2D map of the shock velocity in km s"^ with the contour 
map of the photoionization cone ([0III]/Hy5) overlay ed. Spatial 
scale is 1 arcsec px"^. 



line profiles, 500-800 kms"^ Interestingly, patches of highly 
ionized clumps were observed also outside the ENLR of NGC 
3393 (Cooke et al. 2000) where merging is confirmed by the bi- 
nary black hole system observed in the nucleus (Fabbiano et al. 
2011). Tables [Hand [2] show that where the highest are found 
(e.g. in position 108), the contribution of spectra calculated with 
relatively low velocity - most likely emitted from the ENLR or 
even from the ISIM of the host galaxy - cannot be excluded. The 
relative contribution of the two spectra to each of the observed 
one is revealed by the complex FWHIM of the line profiles cor- 
responding to the diff'erent ions, in agreement with the multi- 
ple kinematics found by Cracco et al. (1201 ll) . For instance, the 
[SII]6716-h lines have a strong contribution from low velocity 
gas. 
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6.3. The distribution of metallicities 



We have calculated the relative abundances to H of the heavy el- 
ements throughout the galaxy, particularly those identified with 
metallicity (e.g. oxygen in the present case), because metallic- 
ity is sensitive not only to star formation, but is also aff'ected 
by the presence of infalls and outflows, i.e. by the interac- 
tion between the merging galaxy and the intergalactic medium 
(Sommariva et al. 201 1). By adopting models which account for 
both shock and photoionization, we found that 0/H is solar in 
almost all positions corresponding to "quiet" gas, if we employ 
the solar relative abundance, 0/H=6.6 x 10""^ reported by Allen 
(1976). However, 0/H becomes about 0.77 solar, when we refer 
to the solar abundance 0/H=8.5 x 10""^ calculated by Anders & 
Grevesse (1989) and Feldman (1992). Sub-solar 0/H values as 
low as 2.6 X 10"^ (Fig. [8]) were found only in some positions far 
from the centre. On the other hand, Cracco et al. (12011b found 
sub-solar 0/H throughout all observed positions using pure pho- 
toionization models. In fact, the SUM A modelling is based on 
both the [OIII] and [Oil] lines. The [OIII] line intensity de- 
pends strongly on the flux from the active center, while the [O II] 
line is afl'ected also by the shock velocity. Moreover, we have 
improved the modelling procedure, accounting for the complete 
spectrum in each position. 

In agreement with Cracco et al. (1201 ID we have found that 
nitrogen is depleted by factors lower than 2 in large regions of 
the cone corresponding to relatively low shock velocities. In the 
regions characterised by high shock velocities (Fig. [5]) N/H is 
roughly solar. S/H shows some depletion in the same external 
south-east regions. N/H is generally high in AGN (e.g. Contini 
et al. l2002a[ Ciroi et al.f2003^. A low N/H may indicate inclusion 
of N into molecules in the ISM and/or acquisition of external gas 
in the merging process. In Sect. [5] we have found that the dust- to- 
gas ratio is relatively low in NGC 7212, showing that dust grains 
and molecules can be destroyed by the itense radiation from the 
AGN and by relative strong shocks. We conclude that sub- solar 
N/H is a further indication of external gas acquisition. 

It is generally known that bursts of star formation are trig- 
gered by galaxy-galaxy mergers (Somerville et al. 2001). As 
mentioned in Sect.[T] Munoz Marm et al. (12007b claim that inter- 
acting galaxies can show UV-light coming from knots produced 
by star clusters. They identify some bright clumps at the south 
of NGC 7212 as star-forming regions. However, the contribution 
of the starbursts to the [O III] 5007 emission line flux is less than 
20%, lower than that of the AGN, in agreement with the results 
for NGC 3393 (Contini 2012), another Seyfert galaxy belonging 
to a merger system. The evidence for star formation motivates us 
to compare the N/0 relative abundance as a function of metal- 
licity calculated throughout the galaxy with the N/0 vs. 0/H di- 
agram presented by Mouhcine & Contini (120021) in Fig. [141 The 
calculated N/0 and 0/H ratios lie in the zone covered by UV 
selected galaxies and massive starbursts in the diagnostic dia- 
gram of Mouhcine & Contini ( 2002, fig.l, black squares and cir- 
cles, respectively). Figure [141 shows a strong variation of N/0 
throughout NGC 7212, indicating a time delay between the re- 
lease of oxygen and that of nitrogen (Contini et al. 2002b ) in the 
regions characterised by 0/H solar abundance. Positions show- 
ing an N/0 overabundance (log(N/0) > -0.86) (see also Fig. [HI) 
may indicate mixing with external matter as a consequence of a 
minor merger. 



S 




-2 - 



7.5 8 8.5 9 9.5 

12 + log(0/H) 

Fig. 14. N/0 vs. 0/H diagram adapted from Mouhcine & Contini 
(2002), (see their fig. 1). The blue area represents a sample of 
UV- selected galaxies, the red area includes starburst nucleus 
galaxies (SBNGs) selected in the optical or in the infrared, and 
the magenta area represents Hii galaxies. The results of our 
model calculations are overplotted (black circles). 



7. Conclusions 

In the previous sections we presented an analysis of the line 
spectra observed by Cracco et al. (2011) in the diff'erent posi- 
tions throughout the ENLR of NGC 7212, adopting calculation 
models which account for the coupled eff'ect of shocks and pho- 
toionization. The results yield new reliable estimates about the 
intensity and distribution of the flux from the AGN, the velocity 
field, the preshock density, the relative abundances, etc. 

In Seyfert galaxies the radiation from the active nucleus is 
the main photoionization mechanism. The AGN in NGC 7212 
is characterised by a flux intensity lower than 10^^ photons 
cm"^ s"^ eV"\ at the Lyman limit, similar to that of low lu- 
minosity AGN (e. g. Cont ini 120041) . rather than to Seyfert 2 
(e.g. Contini et al. I2002al Ciroi et al. 2003). The AGN flux 
peaks in correspondence with the maximum starlight flux. In 
the ENLR of NGC 7212, as in most Seyfert galaxies, the sin- 
gle lines show complex profiles. This indicates that the hydro- 
dynamical field is complex. Shock waves created by winds and 
jets from the stars and/or by collision of matter are generally 
present. Characteristic of shocked clouds is the profile of the 
density downstream, where the gas cools due to recombina- 
tion processes. Compression speeds up recombination, enhanc- 
ing the intensity of lines corresponding to low ionization lev- 
els. The spectra emitted by shocked gas are therefore diff'erent 
from those emitted from clouds heated and ionized by radiation 
from an external source. The relative importance of photoioniza- 
tion and shocks in NGC 7212 throughout most of the 256 spa- 
tial elements simultaneously observed by MPFS, results from 
the consistent modelling of line and continuum spectra. Each 
spectrum stems from a complex grid of models which account 
for the parameter in various ranges. The parameters are con- 
strained by the Hell/H/? line ratios (F), by [S II]6716/6731 (^o), 
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by [O III]: [Oil]: [01] (Vs and F), etc. The shock velocity is also 
constrained by the continuum SED. 

In conclusion, the interaction of galaxy-galaxy or of a galaxy 
with matter from the interstellar or intracluster medium can re- 
veal shocks by line ratios which cannot be readily explained by 
pure photoionization models (e.g. Villar-Martm et al. 1999), by 
high electron temperatures leading to X-ray emission and by ra- 
dio synchrotron radiation. Radiation and collisional processes 
are strongly intermingled in ionizing the gas (e.g. Contini 1997). 
The degree to which shocks contribute to the gas ionization in 
the different regions of NGC 7212 can be found only by consis- 
tent modelling of line and continuum spectra. Our results agree 
with Tadhunter et al. (2000), who claim that shocks contribute 
significantly to the ionization of the ENLR of Seyfert galax- 
ies, even if photoionization contributes substantially in most of 
AGNs (Robinson et al. 1987, Villar-Martm et al. 1997). 
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Table 2. Comparison between calculated and observed line ratios. 



[Oil] [Nelll] Hell [ArlV] [OIII] He I [FeVII] [OI] [Nil] [SII] [SII] 

3727+ 3869+ 4686 4713 5007+ 5876 6087 6300+ 6583+ 6718 6731 



n 


Obs 


Mod 


Obs 


Mod 


Obs 


Mod 


Obs 


Mod 


Obs Mod 


Obs 
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Obs 


Mod 


Obs 


Mod 


Obs 


Mod 


Obs 


Mod 
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35 


2.2 


2.6 


1.23 


1.7 


0.46 


0.58 


0.0 


0.07 


7.43 7.44 
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0.11 


0.0 
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0.45 


4.24 


3.3 


1.08 


0.5 


1.0 


0.6 
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5.2 


4.6 


2.26 


1.1 


0.0 


0.3 


0.0 
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3.57 
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37 
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3. 
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14.6 15. 
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Notes. The + indicates that the doublet is summed up. He II 4686 is blended with [FeIII]4658,4702. 
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Contini et al.: Depletion of heavy elements in NGC 7212 



[Oil] [Nelll] Hell [ArlV] [OIII] He I [FeVII] [01] [Nil] [SII] [SII] 

3727+ 3869+ 4686 4713 5007+ 5876 6087 6300+ 6583+ 6718 6731 
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